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A complete analytical description of the dynamics of current and field driven transverse domain 
walls under the influence of the Dzyaloshinsky-Moriya interaction using the (f) — q model will be 
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parallel or perpendicular to the easy axis anisotropy of the system and a direct reversal respectively 
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I. INTRODUCTION 

Due to the usability for data storage Q, 01 , information 
transport 0] and as logic elements 0,13 magnetic wires 
and spin chains became more and more important during 
the last years. An important role in all these applications 
play the motion of domain walls which underlies some 
criteria: it shall be fast, controllable and stable. 

While the domain wall motion of a normal trans¬ 
verse domain wall is well understood the influence of 
the Dzyaloshinsky-Moriya interaction is mostly tuned 
out. However, the Dzyaloshinsky-Moriya interaction is 
induced by the spin-orbit interaction and appears if the 
symmetry of the system is broken. This is the case at 
surfaces. Now, nearly all magnetic nanowires and spin 
chains are deposit on a surface. Therefore, the influ¬ 
ence of a Dzyaloshinsky-Moriya interaction on the do¬ 
main wall motion should not be ignored. 

A first overview of the influence of the Dzyaloshinsky- 
Moriya interaction on the motion a transverse domain 
wall has been given in the pioneering works of Tretiakov 
and Abanov [g, and Thiaville et al. Q- However, these 
publications don’t cover all possible scenarios. While @ 
describes the influence of a Dzyaloshinsky-Moria vector 
parallel to the easy axis directions in the case of a cur¬ 
rent driven domain wall, 0] describes the influence of 
a Dzyaloshinsky-Moria vector perpendicular to the easy 
axis driven by an external field. Furthermore, in Q the 
description is restricted to the direct reversal scenario 
during the domain wall motion. The possibility of a pre¬ 
cession during the domain wall motion has been ignored. 
Therefore, we can say so far there is no publication which 
describes the field and current driven domain wall mo¬ 
tion under the influence of the Dzyaloshinsky-Moriya in¬ 
teraction, which covers all possible scenarios. This leak 
of information shall be closed with this publication. Fur¬ 
thermore, the description given by Tretiakov and Abanov 
is quite complex. This publication presents a simple way 
which is easy to understand and which offers a deeper 
understanding. 

The next section introduces the q — (j) model which 
has been used to investigate the domain wall motion 


of a transverse domain wall under the influence of the 
Dzyaloshinsky-Moriya interaction. 

II. q — ij> MODEL 

A starting point for the analytical description of any 
kind of domain wall structures and their dynamics is the 
Micromagnetic continuum approximation of the classical 
Heisenberg model. The magnetic properties of a ID or 
quasi-ID domain wall with influence coming from the 
Dzyaloshinsky-Moriya interaction are well described by 
the following Heisenberg Hamiltonian: 

iU) (b> * 

+ ( 1 ) 

i i 

where Si = fl/fis are 3D magnetic moments of unit 
length. 

The first term of this Hamiltonian describes the ferro¬ 
magnetic exchange coupling between nearest neighbors 
with the coupling constant J > 0. The second term de¬ 
scribes the asymmetric exchange (Dzyaloshinsky-Moriya 
interaction) between nearest neighbors coming from the 
spin-orbit coupling. Depending on the symmetry of the 
system the Dzyaloshinsky-Moriya vector Ddm is oriented 
perpendicular or parallel to the bonding axis of the neigh¬ 
boring spins Si and Sj Q- The third sum describes the 
influence of an external magnetic field in z direction 
and the last two terms represent uniaxial anisotropies 
with a hard axis oriented in ±y direction and an easy axis 
oriented in icc direction. Within this manuscript we will 
assume orientations of the Dzyaloshinsky-Moriya vector 
Ddm parallel respectively perpendicular to the easy axis 
anisotropy. The precise orientation will be announced at 
the beginning of every subsection. 

The Hamiltonian Eq. © can be used to calculate the 
dynamics of domain wall numerical (MU. However, a 
discrete structure is not convenient for an analytical de¬ 
scription. In these cases it is comfortable to use the mi- 
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cromagnetic continuum description. The corresponding 
energy: 


+ 00 

E = J £{9,(l))dz (2) 

— OO 

can be got from Ti by p erforming a Taylor expansion 
up to the first order [^. Furthermore, it makes sense 
to use spherical instead of the cartesian coordinates: 
Sx = sin 0 cos 0, Sy = sin0sin(/), and Sz = cosO. The 
detailed energy densities £ corresponding to the different 
scenarios will be given in Sec. IIIII and ITVl 



FIG. 1: Underlying coordinate system: Ke and Kh are the 
uniaxial easy respectively hard axis anisotropy of the system. 
(j> is the angle between the magnetic moment S in the center 
of the domain wall and the x axis; 9 with respect to the z axis 
(not show). The dots symbolizes the next lattice sites of the 
ID or quasi-ID spin system. 

Our goal is to give a complete description of the dy¬ 
namics of a 180° transverse domain wall in a ID or quasi- 
ID system with long axis in z direction. Due to this fact 
we assume for all scenarios that the easy axis anisotropy 
is oriented in ±z direction to take into account the con¬ 
tributions from crystalline anisotropy as well as dipole- 
dipole interaction [9|. Furthermore, if Dh ^ 0 (additional 
hard axis anisotropy) the hard axis is oriented in ±?/ di¬ 
rection (see Fig.[T]). Under these assumptions in all cases 
the magnetic moments in the domains are oriented in ±z 
direction and the orientation of the magnetic moment in 
the center of the domain wall depends on Dh'. if Dh > 0 
the magnetic moment will be oriented in one of the x 
directions (ix, (/) = 0 or (/) = tt), if Dh = 0 the symmetry 
in the xy-plane is not broken and the magnetic moment 
can be oriented in any of the possible directions in the 
xy-plane (0 < (/) < 27r). So far we haven’t take into 
account the Dzyaloshinsky-Moriya interaction. In the 
following we assume that the Dzyaloshinsky-Moriya in¬ 
teraction is small with respect to the other energy contri¬ 
butions (exchange and anisotropy): |Ddm| = Ddm ^ J, 


Ddm De, and Ddm ‘C Dh- The assumption of a 
small Dzyaloshinsky-Moriya interaction is necessary to 
find domain wall solutions instead of spin spirals. 

Within the micromagnetic continuum description it is 
easy to verify that with the assumptions made before the 
static solution (Bz = 0) of a 180° transverse domain wall 
is given by |12 h14| : 

( z — q{t) \ 

e j , (3) 

and 

Sz = COS0 = ±tanh 

This profile is the general solution of a ID or quasi- 
ID 180° domain wall and independent from the explicit 
form of the Hamiltonian given by Eq. du, as long as the 
Dzyaloshinsky-Moriya interaction is small: 0 is the an¬ 
gle with respect to the z axis, A the domain wall width, 
which depends on "H, and q{t) is the center of the domain 
wall. The center of the domain wall itself is character¬ 
ized by Sz = cos 0 = 0 which means 0 = 1:12. Additional, 
we find for the variation b£ 150 of the microscopic energy 
density £ under the condition to be in the center of the 
domain wall {0 = 7r/2): 

= MsBz . ( 5 ) 

e=7r/2 

A detailed description is not necessary at this point and 
will be given later in Sec. IIIII and IIVI As before in Eq. ([T]) 
Bz is the external field. Mg = is the saturation 

magnetization, and a is the lattice constant. Additional 
we can deduce the following conditions from the domain 
wall profile: 


J0 


d0 sin0 , , 

dl = <•=> 

and 

dt sin0 dt ’ 

where v is the velocity of the domain wall. 

The underlying equation of motion is the Gilbert equa¬ 
tion with additional spin torque terms (adiabatic and 
non-adiabatic) describing the influence of an electric cur¬ 
rent Uz oc J, with j the current density Eini: 


dt 


1 Q 

7^ X Heff - X — 

dt 

dS ^ ^ dS 

Uz- -h pUzS X 


dz 


dz 


( 8 ) 


The first term of this equation describe the precessional 
motion of the magnetic moment S in the effective field 
Heff = -{i/ps){dn/dS). 7 is the gyromagnetic ratio. 
The second term is the so called Gilbert damping with 
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damping constant a [l5l| . The third and fourth terms 
are the adiabatic and non-adiabatic spin torque terms 
describing the influence of an polarized electric current. 
The adiabatic spin torque term appears due to the mo¬ 
mentum conservation during the interaction between the 
localized magnetic moments S and the spins of the elec¬ 
trons. The non-adiabatic spin torque term takes into ac¬ 
count the scattering processes of the electrons. /3 is the so 
called non-adiabaticity constant. For a detailed descrip¬ 
tion see e.g. The corresponding micromagnetic 

equations in spherical coordinates are: 


dt 


7 . dcj) de . d(j) 

77—^77 -asm6»— -u^- - /du^sinO— 

MssmOocj) dt dz dz 


(9) 


and 

d(j) 5£ a dO dcf) /3uz d9 

dt Ms sin 0 60 sin 9 dt ^ dz sin 9 dz 

( 10 ) 


S£/69 respectively 6£ / 6(j) are the variations of the energy 
density with respect to 9 respectively (j). 

Inserting the conditions Eq. together with the 

assumption to be in the center of the domain wall: 
9 = 7r/2, lead to the following equations which are the 
starting point of our considerations: 


V 

A 


j 6£ dcj) Uz . d^ 

Ms 6(1} “ dt A dz ’ 


( 11 ) 


and 


d(j) 

dt 


n aw 

-jBz + -^ -Uz 


d(j) 

dz 


I3uz 

A 


( 12 ) 


The first equation m describes the velocity v and the 
second equation m the precession of the domain wall. 

These equations are simple to solve and depend only on 
q (or V = dq/dt) and (j). In the next sections the explicit 
form of (j}{z, t) under the assumptions of a Dzyaloshinsky- 
Moriya interaction with a vector Udm parallel and per¬ 
pendicular to the easy axis anisotropy shall be charac¬ 
terized and the corresponding equations (EH) and EH 
solved. 

A detailed outline of the following sections can be 
found in following table: the first and second column give 
the section and sketch (Fig.) of one scenario and the 
last four column the corresponding informations about 
the used assumptions (orientation of Udm, Dg and Dh 
as well as the reversal mechanism where direct means 
d(j)/dt = 0 and precession d(j)/dt ^ 0). 

This table can be seen as a road map of this publi¬ 
cation. The calculations themselves (dealing with main 
equations Eq. EH and EH) are simple and follow al¬ 
ways the same procedure: After writing down the energy 


TABLE I: 


Sec. 

Fig. 

De 

Dh 

Ddm 

reversal 

Ill A 

|2| 



+3 

direct 

IIIB 

|3| 


no anisotropy 

+z 

precession 

JVAI 

m 


±y 

+y 

direct 


p 



+X 

direct 

^vcl 


±2 

no anisotropy 

+y 

precession 


density £, corresponding to % [Eq. (P)] and the made 
assumptions (see table |I|, the variation 6 £/ 6 (j)\e=-K /2 = 0 
has to be calculated. Inserting the result in the main 
equations Eq. EH and EH together with an assumption 
about the angle (p lead to the hnal equations describ¬ 
ing the domain wall motion. Independent of the chosen 
assumptions 6£I69 = 0 lead to the domain wall profile 
Eq. (|4]) together with the domain wall width A. 


III. DM VECTOR PARALLEL TO THE EASY 
AXIS ANISOTROPY DIRECTIONS 


In the following we assume a spin chain with a trans¬ 
verse domain wall with an easy axis anisotropy in dzz 
direction. We will stay with this assumption till the end 
of this publication. Furthermore, in this section we as¬ 
sume that the Dzyaloshinsky-Moriya vector is parallel to 
the easy axis direction: Ddm = Dy)m£z- Then in the fol¬ 
lowing two subsections the two possible scenarios Dh ^ 0 
and Dh — 0 will be discussed. In the first scenario the 
hard axis anisotropy Dh leads to a symmetry break and 
a direct reversal of the magnetic moments during the 
motion. In the second scenario (no hard axis anisotropy, 
Dh = 0) the system shows a rotational symmetry around 
the long axis {z axis). 

In both cases the micromagnetic energy density £ is 
given by: 


£ = A 


d9 


dz 


— + sm" 9 — 


d(j) 


dz 




-I- Ddm sim 9— 
dz 


— MsBz cos9 + KhSrn^ 0 sin^ (p — KeCOS^ 9 , (13) 


where A = J/2a is the exchange stiffness, AT,, = 

(j] G {e,h}) the uniaxial anisotropies, Bz the external 
field, Vdm = Ddm/o? the micromagnetic Dzyaloshinsky- 
Moriya interaction, Ms = t^s/o? the magnetization, and 
a the lattice constant. The difference between both sce¬ 
narios is the appearance of the hard axis anisotropy: first 
scenario Kh ^ 0 (direct reversal), second scenario Kh = 0 
(processional motion). 

In the following we will continue with the micromag¬ 
netic description instead of the discrete description used 
in Eq. (P. 
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A. Direct Reversal 

In this subsection we assume a hard axis anisotropy 
Kfi 7 ^ 0 in ±y direction which is dominating: ^ 

This immediately leads to the fact that we can 
expect: di^/dz = 0. In principle due to Kh T’dm the 
Dzyaloshinsky-Moriya interaction has no influence and 
we can expect in this case the well known results of do¬ 
main wall motion without Dzyaloshinsky-Moriya interac¬ 
tion. 



FIG. 2: Sketch of the domain wall profile of a ID trans¬ 
verse domain wall with Dzyaloshinsky-Moriya interaction cor¬ 
responding to the scenario described in Sec. IIII Al left hand 
side X — z plane, right hand side (circle) x — y plane. The 
Dzyaloshinsky-Moriya vector 2?dm is oriented parallel to the 
easy axis and perpendicular to the hard axis anisotropy 
Kh. Due to the hard axis anisotropy Kh we can assume a 
direct reversal A(f)/At = 0. 


It is quite easy to verify that in this case 5EI59 = 0 
if we assume Bz = 0 (static solution) for this moment 
leads to the well known domain wall profile (| 31 ) , with the 
domain wall width: 


A = 


A 


Ke + Kh sin^ 4> 


(14) 


Then, with under the assumption to be in the 

center of the domain wall: cos0 = 0 => d = 7 r/ 2 : 




Khsu\{2(j)) , 


the equations CD and m become: 


V 

A 


Khl 

Ms 


sin(2())) -I- 


Uz 

'a 


and 


av 

0 = -jBz + 


PUz 

A 


(15) 


(16) 


(17) 


Here, the domain wall width A is given by Eq. (fHl) . 

Inserting Eq. CD in (iroi) leads to the condition which 
describes the stability of the domain wall (Walker break¬ 
down) : 


M.s ’ 


= —arcsin 
2 


MsBz + {P-a )^l^ 
aKh 


(18) 


Inserting this result in Eq. CD, gives the domain wall 
width A during the domain wall motion e.g. Bz ^ 0. 
Furthermore, Eq. dT51) together with Eq. ((T5)) leads to the 


well known velocity e quat ion for a domain wall motion 
with direct reversal IllL [l8l| : 

V =-A -I- Uz . (19) 

a a 

The first term describes the influence of an external field 
and the second term the influence of the electric cur¬ 
rent. The same result can be derived also directly from 

Eq. (fT71) . 

As discussed before the dominating hard axis 
anisotropy Kh leads to the fact that the Dzyaloshinsky- 
Moriya interaction has no influence. This can be seen 
also in the equation for the velocity Eq. CD and sta¬ 
bility equation (fT51) . The equations for the velocity as 
well as the stability of the domain wall are the same as 
in this scenario. For comparison, the domain wall mo¬ 
tion of a transverse domain wall without influence of a 
Dzyaloshinsky-Moriya interaction is described in [Hill. 

On the other hand the assumption of a dominating 
hard axis anisotropy Kh is a quite hard criteria which 
was necessary to solve this problem analytical. A softer 
assumption: Kh not dominating will lead to an influence 
of the Dzyaloshinsky-Moriya interaction. However, 
in this case a complete analytical description without 
approximations is impossible, e.g. Tretiakov and Abanov 
0 have used a perturbation theory to obtain their result. 


B. Precessional motion 

While Kh 7 ^ 0 is not exactly solvable without the 
assumption of a dominating hard axis anisotropy the 
situation changes if we assume Kh = 0 (no hard axis 
anisotropy). As before in previous subsection we assume 
the Dzyaloshinsky-Moriya vector parallel to the easy axis 
anisotropy in -\-z direction. Furthermore, we assume 
Kh = 0 , therefore we have a rotational symmetry around 
the z axis and this means we cannot neglect the influence 
of the Dzyaloshinsky-Moriya interaction: d(j)/dz 7 ^ 0. 
Furthermore, we have to expect that the magnetic mo¬ 
ments and therefore the domain wall process during the 
motion: dp/dt 7 ^ 0 . 



FIG. 3: Sketch of the domain wall profile of a ID trans¬ 
verse domain wall with Dzyaloshinsky-Moriya interaction cor¬ 
responding to the scenario described in Sec. IIII Bl left hand 
side X — z plane, right hand side (circle) x — y plane. The 
Dzyaloshinsky-Moriya vector 22 dm is oriented parallel to the 
easy axis anisotropy K^. Dne to the Dzyaloshinsky-Moriya in¬ 
teraction and the absence of the hard axis anisotropy Kh = 0 
the domain wall profile shows a twist and we have to assume 
a precessional motion A4>/At A 0. 
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In this scenario 5E jdcj) = 0 leads to: 


54) 


_d 

'dz 


sin^ 6 ( 2A-^ ± 2?dm 

dz 


= 0 , 


( 20 ) 


which means: 


X>dm 
^ 2A 


( 21 ) 


The solution of this differential equation is given by: 


1) = Tr(z - qit)) + 4)0 , (22) 


where we have set zq = q and T = 2 ?dm/ 2^- This 
equation describes the linear increase / decrease of ^ 
with 0 due to the Dzyaloshinsky-Moriya interaction. In 
other words we find a twist in the domain wall profile: 
d4>/dz ^ 0 (see Fig. [3]). 

With this knowledge it is easy to show that SSjSO = 0, 
with Bz = d, leads to the well known domain wall profile 
(|3|) together with the domain wall width: 


A = 



(23) 


To describe the dynamics it is necessary to know 5E/d4) 
and d4>/dz: SE/S4> is equal to zero [see Eq. (1301) ] and 
d4)/dz is described by Eq. (1311) . With this information 
we can write down our main equations Eq. m and © 
as: 


are identical. However, due to the twisted shape of the 
domain wall not only the domain wall rotates during the 
motion but also our coordinate system. This rotation 
comes from the twisted shape but does not take into ac¬ 
count the precessional motion of the magnetic moments. 
In other words Eq. (1371) leads to wrong results because 
we are not in a stationary frame or in the rotating frame 
with the correct rotational frequency. In the later case we 
can expect d4>/dt = 0. So, if we are interested in investi¬ 
gating the oscillation of the domain wall we have to be in 
a stationary or rotating frame with the same rotational 
frequency as the magnetic moments. However, both is 
not the case. Therefore, we have to make an additional 
transformation to get in one of these two frames. In the 
following we use the stationary frame. 

We know from Eq. ((3T|) that: 

^ = =Fr O d(/) = Thdz , (28) 


and therefore: 


# ^ Y— 

dt dt 


(29) 


This is the correction we have to take into account 
to change from the rotating to the stationary frame. 
V = dz/dt is the velocity of the domain wall and at the 
same time the velocity of the moving frame. Adding this 
term to Eq. (1371) and taking into account Eq. (133|) . the 
precession in the stationary frame is described by: 


V _l_ « T- 


d4) 
dt 

Inserting (1351) in 


av rj , r 


we get: 


(24) 

(25) 


iBz l + aP ^ a- P 

V = -+ ■;-T ^-?rl Auz . 

V 1 + a 


cx 


1 + a^ 


(26) 


This result is identical to the velocity of a field and cur¬ 
rent driven domain wall without Dzyaloshinsky-Moriya 
interaction [ll|, [l^, just with the additional term 
=F(q! — /3)/(l + a‘^)TAuz describing the influence of the 
Dzyaloshinsky-Moriya interaction. This term increases / 
decreases the velocity depending on the direction of the 
Dzyaloshinsky-Moriya vector. 

Inserting in (1351) leads to: 


d0 a — f3uz 'jBz . (l-l-a/3)r 

_ —_— _ ^ _ 'll 

dt 1 -I- A 1 -I- 1 -I- 


This equation describes the precession of the transverse 
component of the domain wall during the motion. 

At this point it should be noticed that we have assumed 
to be in the center of the domain wall. This means that 
we are in a moving frame. In the case of the domain 
wall velocity this is not a problem, because the velocity 
of the domain wall and the velocity of the moving frame 


d4) 

dt 


a — P 
1 -I- 



iT^A 


iBz 
1 -I- 


(liaTA) . 


(30) 

Together with equation (1351) this equation describes the 
precessional domain wall motion of transverse domain 
walls with a Dzyaloshinsky-Moriya interaction where the 
Dzyaloshinsky-Moriya vector is parallel to the easy axis 
anisotropy. This result is identical with the result given 
by Tretiakov and Abanov Q , however the given descrip¬ 
tion here is simpler. 


IV. DZYALOSHINSKY-MORIYA VECTOR 
PERPENDICULAR TO THE EASY AXIS 
ANISOTROPY DIRECTIONS 

The assumption of the previous section was that the 
Dzyaloshinsky-Moriya vector is parallel to the easy axis 
anisotropy. In principle this not the general situation. 
Theoretically, the Dzyaloshinsky-Moriya vector can show 
in any direction, depending on the symmetry of the sys¬ 
tem. To simplify the problem and due to the fact that 
we are interested in analytical solvable problems we re¬ 
strict ourself to scenarios with the Dzyaloshinsky-Moriya 
vector perpendicular to the easy axis anisotropy. In 
the following, we investigate three scenarios: (1) the 
Dzyaloshinsky-Moriya vector has an orientation perpen¬ 
dicular to the easy axis anisotropy K,, and parallel to 
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the hard axis anisotropy Kh, (2) the Dzyaloshinsky- 
Moriya vector is perpendicular to the easy as well hard 
axis anisotropy Kh-, and (3) there is just an easy axis 
anisotropy {Kh = 0). In all these scenarios we assume 
that the spin chain has an alignment along the 0 axis, the 
easy axis is oriented in ±z direction (transverse domain 
wall) and the Dzyaloshinsky-Moriya vector is oriented in 
y direction. To distinguish the three scenarios we change 
the strength (1),(2) Kh ^ 0, (3) Kh = 0 and orienta¬ 
tion of the hard axis anisotropy Kh'. (1) the hard axis 
anisotropy Kh is oriented in ±y direction, (2) in ±x di¬ 
rection. 

The micromagnetic energy density £ in the case of a 
Dzyaloshinsky-Moriya vector in y direction parallel to the 
hard axis anisotropy Kh is given by: 


£ = A 




X>dm cos ( j ) 


di 


— MsB^cosO + KhSin^ (f> — KgOOS^ 9 , (31) 


For a hard axis anisotropy Kh in ±x direction we have 
to replace the term Kh sin^ 6 sin^ (j) by Kh sin^ 9 cos^ (j). 


A. Direct Reversal: Dzyaloshinsky-Moriya vector 
parallel to Kh 


As before in the first scenario, we assume the appear¬ 
ance of both uniaxial anisotropies Ke and Kh- However, 
this time we assume that the Dzyaloshinsky-Moriya vec¬ 
tor Vdm shows in +y direction (parallel to Kh) and not 
in +x direction. In this case we don’t need to assume 
a dominating hard axis anisotropy to solve the problem 
exactly. 

In this scenario the variation 6£/69, with = 0, leads 
to the following result: 


J9 


2A-— - X>dm sin ( j )— 
dz^ dz 


2 sin 9 cos 9 


Ke + Kh sin^ 6 — A 



(32) 

= 0 



Ke 


© © 
■Udm Kk 



FIG. 4: Sketch of the domain wall profile of a ID trans¬ 
verse domain wall with Dzyaloshinsky-Moriya interaction cor¬ 
responding to the scenario described in Sec. IIV Al left hand 
side X — z plane, right hand side (circle) x — y plane. The 
Dzyaloshinsky-Moriya vector Odm is oriented perpendicular 
to the easy axis Ke and parallel to the hard axis anisotropy 
Kh- The Dzyaloshinsky-Moriya interaction only affects the 
domain wall width but not the profile itself. Due to the hard 
axis anisotropy Kh we can assume a direct reversal d4>/dt = 0. 


Therefore, we are able to write the main equations 


Eq. ITT]) and 

m 

as: 


V 

A 

7 

Ms 

1 2Kh cosp+ ^ j sm P + ^ 

(35) 

0 = 

av 

~A ■ 


(36) 

After eliminating 

V we get the stability condition: 




^Bz + iP-a)^ 

(37) 

sm(p ■■ 

^{2KhCOsP+^) ■ 


Inserting this result in (1551) gives the formula for the ve¬ 
locity: 

iBz * P 

V = -A H- Uz ■ (38) 

a a 

This result is identical to the velocity of a transverse 
domain wall with direct reversal and no Dzyaloshinsky- 
Moriya interaction [l^, [l^ . The only difference between 
both descriptions is the stability criteria (|37)) . In the case 
with Dzyaloshinsky-Moriya interaction the domain walls 
are more stable, which means the Walker breakdown 
appears at higher field and current values. 


B. Direct Reversal: Dzyaloshinsky-Moriya vector 
perpendicular to Kh 


With the assumption that (p is constant: dp/dz = 0 and 
d(j)/dt = 0 this differential equation can be solved easily 
and leads again to the well known domain wall profile ([3]) 
together with an domain wall width which is not influ¬ 
enced by the Dzyaloshinsky-Moriya interaction: 


A = 


A 


Ke + Kh sin^ (p 


(33) 


The second variation 5 £/under the assumption to be 
in the center of the domain wall: 9 = tt 12 and a constant 
p, leads to: 

^ = (^2KhCOSp+'^^^ sinp . (34) 


In the following we assume the same geometry as 
before, however with the following difference: the 
Dzyaloshinsky-Moriya vector is now oriented in +x di¬ 
rection and therefore perpendicular to both anisotropies 
Ke and Kh (see Fig. [5]). 

This scenario is solvable if we assume that 
both anisotropies are dominating with respect to 
Dzyaloshinsky-Moriya interaction: Kh S> ©’dm, Ke S> 
Ddm • In this case the Dzyaloshinsky-Moriya interaction 
can be neglected and we find a normal transverse domain 
wall with the profile ([3]) and the domain wall width: 


A = 


A 


Ke + Kh sin^ p 


(39) 
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FIG. 5: Sketch of the domain wall profile of a ID trans¬ 
verse domain wall with Dzyaloshinsky-Moriya interaction cor¬ 
responding to the scenario described in Sec. IIV Bl left hand 
side X — z plane, right hand side (circle) x — y plane. The 
Dzyaloshinsky-Moriya vector 2?dm is oriented perpendicu¬ 
lar to the easy axis and hard axis anisotropy Kh- The 
Dzyaloshinsky-Moriya interaction shows no influence on the 
profile in this case. Due to the hard axis anisotropy Kh we 
can assume a direct reversal A(f>/At = 0 . 


Then, the stability is characterized by: 


1 

= -arcsin 


aKh 


(40) 


and the velocity of the domain wall described by: 


iBz , P 

V = -A -I- Uz 

a a 


(41) 


C. Precessional motion 

In this subsection we assume Kh = 0, which means 
no hard axis anisotropy and therefore a precessional mo¬ 
tion d(j)/dt ^ 0. The Dzyaloshinsky-Moriya vector in 
this scenario is oriented in direction (see Fig. [S]). 
Furthermore, for simplicity we assume that dp/dz = 0. 



FIG. 6 : Sketch of the domain wall profile of a ID trans¬ 
verse domain wall with Dzyaloshinsky-Moriya interaction cor¬ 
responding to the scenario described in Sec. IIV Gl left hand 
side X — z plane, right hand side (circle) x — y plane. The 
Dzyaloshinsky-Moriya vector Ddm is oriented perpendicular 
to the easy axis anisotropy Tfe. Due to the absence of the 
hard axis anisotropy Kh we can assume a precessional mo¬ 
tion d(l>/dt ^ 0 which leads to a time dependent domain wall 
width A(t) and domain wall energy i?(t). 

This is definitively the case when the transversal compo¬ 
nent of the domain wall and the Dzyaloshinsky-Moriya 
vector Ddm are perpendicular as in Fig. [51 In this case 
the Dzyaloshinsky-Moriya interaction assists the magne¬ 
tization reversal in the x — z plane which is described by 
the angle 0. The assumption dp/dz = 0 is not neces¬ 
sarily the case when Ddm is parallel to the transversal 
component of the domain wall. 


Under these assumptions, especially dp/dz = 0, we 
find: 

5E (^0 

— = 2A—^ -I- 2Ke sin 0 cos 0 — 0 , (42) 

00 dz^ 

and therefore, the well know domain wall profile (|3]) to¬ 
gether with the domain wall width: 


A M 

V ’ 

and the domain wall energy: 

E = Ay/AKe + ttHdm ■ 


(43) 


(44) 


The variation S£/6(j) under the assumption to be in the 
center of the domain wall {0 = 7 r/ 2 ) leads to: 


S(j) 


Bum 


■ sm ( 


e=7r/2 


(45) 


Then, it is easy to write the equations m and m as: 

V dp J Bum ■ , Uz . . 

-r = a--—- 7 — sin 0 4 - — , (46) 

A dt Ms A A ^ ^ 

and 


d 0 av Puz , , 

= (47) 

Eliminating dp/dt leads to the following velocity equa¬ 
tion: 


"fBz . , l + aP iBum sin p{t) 
v{t) = -tA -I- —-Uz - 


a + - 


1 


MsA(l-ba2) ■ 


(48) 


In time average (mean value) the velocity of the do¬ 
main wall is equal to: 


iBz . 

V = -r A ■ 


1 + ap 


(49) 


which is equal to the velocity of a transverse domain wall 
with precession and without Dzyaloshinsky-Moriya inter¬ 
action m- 

Eliminating v in Eq. (H51) and (ITTl) lead to the following 
differential equation for p: 


dp 

dt 


a — a Uz 


iBz 


aiBuM 


1 -|- A 1 -|- M 5 A (1 -|- a^) 


sin0(t) . (50) 


This differential equation is needed to be solved to cal¬ 
culate the velocity v(t). For simplification we write the 
differential equation Eq. dSOl) as: 


dp / • , / \ 

— = 00-01 sin 0 (t) , 
dt 


with 


Po — 


a — a Uz 


iBz 


A 1 + 


(51) 


(52) 
























FIG. 7: Domain wall motion of a transverse domain wall with 
Dzyaloshinsky-Moriya interaction corresponding to the sce¬ 
nario described in Sec. II V Cl domain wall position as function 
of time corresponding to the scenario that the Dzyaloshinsky- 
Moriya vector Udm shows an orientation perpendicular to the 
easy axis anisotropy and the assumption that there is no 
additional hard axis anisotropy Kh = 0. The last assumption 
leads to a processional motion which leads to the oscillation 
in z{t). The curves coming from computer simulations with: 
HsBz/J — 0.2, Dom/J = 0.004, and De/J = 0.005. 


the length of 1000 magnetic moments. The orientation 
of the chain is along the z axis and the magnetic momets 
inside the domains are in ±z direction. The unrelaxed 
domain wall is a 180° head-to-head transverse domain 
wall which becomes the final shape during the relaxation 
process. After switching on an electric current the do¬ 
main wall displacement has been determined by looking 
for the zero-crossing of the tanh-profile of the domain 
wall: Sz = 0- For more details of the computer simula¬ 
tions on the atomic length scale please see [§-[ll|. 

Fig.[7]shows the domain wall displacement as function 
of time for a current driven transverse domain wall de¬ 
scribed as in the last scenario with a < fi. As predicted 
by analytical calculations the domain wall shows an os¬ 
cillatory behavior and moves with the averaged velocity 
given by Eq. (|4^ . However due to the oscillation the 
domain wall arrives the end of the system with a delay 
or earlier depending on the sign of the Dzyaloshinsky- 
Moriya vector X>dm ■ 


V. SUMMARY 


and 


4>i 


MsA(l -ba2) ■ 


(53) 


Eq. m can be easily solved by separation of variables. 
The result is: 


(/)(t) = 2arctan 


< 1^1 


tan 




(54) 

An interesting situation appears if the domain wall is 
driven by an electrical current [uz 0,Bz = 0) and if 
we assume a = j3. It is known that in this case a trans¬ 
verse domain wall without Dzyaloshinsky-Moriya inter¬ 
action shows no precession difi/dt = 0 even if the domain 
wall normally process difi/dt ^ 0 if a /3. However, in 
the case with Dzyaloshinsky-Moriya interaction we find a 
processional motion even in this special situation: a = /3, 
Bz = 0. In this case we find: 


(j){t) = 2arctan 


exp 




A/sA(l -ba2) 



(55) 


and 


V = u 


z ■ 


(56) 


To get a better feeling about the accuracy and cor¬ 
rectness of the results of the last subsection additional 
numerical simulations have been performed. Therefore 
the Gilbert equation ([8]) has been solved using the dis¬ 
crete Heisenberg model (HD and as numerical solver the 
Heun method. Starting point of the simulation was a 
relaxed domain wall in the first half of a linear chain of 


In summary this manuscript gives a complete overview 
of the domain wall motion of transverse domain walls 
under the influence of the Dzyaloshinsky-Moriya interac¬ 
tion. The description is analytical using the q — (j) model 
and covers the scenarios of a Dzyaloshinsky-Moriya vec¬ 
tor parallel respectively perpendicular to the easy axis 
direction and the fact that the domain wall can show a 
direct reversal or processional motion. The first two sce¬ 
narios describe the situation of a Dzyaloshinsky-Moriya 
vector parallel to the easy axis direction. These scenarios 
have been described by Tretiakov and Abanov @ how¬ 
ever using too complex description. The given descrip¬ 
tion within this paper is simpler and therefore easier to 
understand. The third scenario is identical to the de¬ 
scription of Thiaville et al. Q of a Dzyaloshinsky-Moriya 
vector perpendicular to easy axis and parallel to the hard 
axis direction. However, here we give some more infor¬ 
mations. The fourth scenario describes the situation of 
a Dzyaloshinsky-Moriya vector perpendicular to the easy 
axis and at the same time hard axis anisotropy. In this 
case the Dzyaloshinsky-Moriya interaction has no influ¬ 
ence as long we assume that X>dm is small with respect to 
the anisotropies Kg and Kh- The fifth scenario assumes a 
Dzyaloshinsky-Moriya vector perpendicular to easy axis 
direction. However, within this scenario we assume a pre- 
cessional motion. This scenario has not been described 
so far and compliments the description of the dynam¬ 
ics of transverse domain walls under the influence of a 
Dzyaloshinsky-Moriya interaction. We have seen that 
the Dzyaloshinsky-Moriya interaction can lead to a mod¬ 
ification of the velocity, or at least the stability condition. 
Other scenarios show no influence of the Dzyaloshinsky- 
Moriya interaction on the velocity. To compliment these 
results additional computer simulations could be help¬ 
ful, especially in cases where the domain wall cannot be 
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described by a simple transverse domain wall. 
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